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ABSTRACT 

Context. Low- and intermediate-mass stars go through a period of intense mass-loss at the end of their lives during the asymptotic 
giant branch (AGB) phase. While on the AGB a significant part, or even most, of their initial mass is expelled in a stellar wind. 
This process controls the final stages of the evolution of these stars and contributes to the chemical evolution of galaxies. However, 
the wind-driving mechanism of AGB stars is not yet well understood, especially so for oxygen-rich sources. Characterizing both the 
present-day mass-loss rate and wind structure and the evolution of the mass-loss rate of such stars is paramount to advancing our 
understanding of this processes. 

Aims. We study the dusty wind of the oxygen-rich AGB star W Hya to understand its composition and structure and shed light on the 
mass-loss mechanism. 

Methods. We modelled the dust envelope of W Hya using an advanced radiative transfer code. We analysed our dust model in the 
light of a previously calculated gas-phase wind model and compared it with measurements available in the literature, such as infrared 
spectra, infrared images, and optical scattered light fractions. 

Results. We find that the dust spectrum of W Hya can partly be explained by a gravitationally bound dust shell that probably is 
responsible for most of the amorphous AI 2 O 3 emission. The composition of the large (~ 0.3 pm) grains needed to explain the scattered 
light cannot be constrained, but probably is dominated by silicates. Silicate emission in the thermal infrared was found to originate 
from beyond 40 AU from the star. In our model, the silicates need to have substantial near-infrared opacities to be visible at such large 
distances. The increase in near-infrared opacity of the dust at these distances roughly coincides with a sudden increase in expansion 
velocity as deduced from the gas-phase CO lines. The dust envelope of W Hya probably contains an important amount of calcium 
but we were not able to obtain a dust model that reproduces the observed emission while respecting the limit set by the gas mass-loss 
rate. Finally, the recent mass loss of W Hya is confirmed to be highly variable and we identify a strong peak in the mass-loss rate that 
occurred about 3500 years ago and lasted for a few hundred years. 


1. Introduction 

Low- and intermediate-mass stars in the last stages of their 
lives evolve to the asymptotic giant branch (AGB) and occupy 
the region of the Hertzsprung-Russell diagram of high lumi¬ 
nosities (typically 5000 - 50000 Lq) and low effective tem¬ 
peratures (~3000K). Dust grains can form in the outer layers 
of the extended atmospheres of AGB stars where the temper¬ 
atures are low and the densities are enhanced by stellar pul¬ 
sations. Radiation pressure acting on the grains causes a slow 
wind to develop lOkm/s) through which the star looses mass, 
a process that controls th e subsequent evolution of the star 
(iHabing & 01ofssonll2003l) . The cumulative mass lost during the 
AGB phase is important for the chemical evolution of galax¬ 
ies, as the atmospheric composition of AGB stars is modified 
by dredge-ups, which brin g elements sy nthesized in the stellar 
interior to the surface (e.g. lHabingll 19961) . 

Despite a qualitative understanding of the AGB mass-loss 
process, models are not yet able to predict the mass-loss history 
of an AGB star from first principles. Both the magnitude of the 
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mass loss and its variation in time are still unanswered ques¬ 
tions. For instance, model calculations show that the winds of 
oxygen-rich AGB stars cannot be driven by absorption of pho¬ 
tons. The reason is that iron-bearing silicates that have large 
near-infrared absorption cross-sections to acquire enough mo¬ 
mentum also heat up ve ry efficiently and cannot exist in the 
wind acceleration reg ion (IW oitkell2006l). A solution to this prob¬ 
lem was proposed bv lHofne" ( 20081) . who suggested that refrac¬ 
tory species with small near-infrared absorption cross-sections 
- that can exist close to the star - can still transfer momentum 
to the wind by means of scattering, if the grains can grow to 
a few tenths of a micron in diameter. Indeed, in a more com¬ 
prehensive study of the possi ble wind-drivers in oxygen-rich 
AGB stars, iBladh et^ (1201 3l) found that the dust species ex¬ 
pected to form in an oxygen-rich environment are not able to 
drive the wind via the absorption of photons. The nature of 
the grains responsible for driving the wind and the particulars 
of the wind-drivin g mechanism are still widely discussed (e.g . 
Norris et alJ I 20 II iBladh & Hofnd |2^ iBladh et al.1 |2^ 
Gail et alJ 1201 3t iKarovicova et alJ 2013 ). From wind-driving 


calcula tions, the preferred wind-driving cand idates are silicate 
grains (iHofned [ZOOSt IBladh & Hofneill2012h . However, amor¬ 
phous aluminum oxide is also found to be an important com- 
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ponent in general oxygen-rich AGB stars with low mass-loss 


rates (e.g. iLorenz-Martins & Pomneia 

20001; LSneck et al.ll200fll 

iHeras & Honvl I 2 OO 5 I: Karovicova et al. 

I 2 OI 3 I). The latter might 


act as a seed for further grain growth and be crucial for start¬ 
ing the wind, despite their low abundance. Aluminum oxide has 
been suggested to form in the extended atmospheres of oxygen- 
rich late-typ e stars both on the basis of grain con densation calcu¬ 
lations (e.g. W)itkdl2QQf ) and o bservations (e.g. lVerhoelst et all 
l2006t IZhao-Geisler et al, 1201 Ih . MIDI visibilities of RT Vir, for 


instance, point t o the presence of am orphous AI 2 O 3 close to the 
star (at ~ 1.5 R.i.: ISacuto et al.ll2013h . However, dust excess spec¬ 
tral fits often require a too large amount of this grain species 
to match the observed infrared emission; so large that it often 
implies a super-solar a luminum abundance in the outflow (e.g. 
iKarovicova et al.ir2013h . 

We model the dust envelope of W Hya, a clos e-by oxygen- 
rich A GB star with a fairly low mass-loss rate (iKhouri et al.l 
l2014aL henceforth Paper I) for which a wealth of observa¬ 
tional data is available. These data probe the wind structure 
of W Hya on spatial scales ranging from tens of mas to min¬ 
utes of arc and reveal grains with radii of about 0.3 iim a t 
40 mas (or 2R*) seen in scattered light (iNorris et al.ll2012h : 
an inner radiu s for the amorphous sili cate emission of 500 
mas (or 25 R*; IZhao-Geisler et all 1201 ih . and a mass loss that 
is variable on d ynamical timesc ales of thousands of years 
dHawkinsI Il990t fCox et af] l 2012 h that resembles the forma¬ 
tion of detache d-shells (exclusively) seen around carbon-rich 
AGB stars (e.g. lOlofsson et aT]ll99(]l [l996f) . Although WHya 
has previously been studied extensively, the models obtained 
for the dust exc ess only focused on t he emission from the 
inner wind (e.g. Heras & Honvl l2005t lJusttanont et al.l l2004t 
iKhouri et alJ2014ah. None ha ve accounted for the scattered light 
fractions seen by Norris et^ or t he inner emission ra dius of the 
amorphous silicates reported by IZhao-Geisler et all however. 
Furthermore, all dust emission models obtained so far require 
super-solar aluminum abundances in the dust when compared to 
the gas-phase mass-loss rate. 


Based on the strong evidence for amorphous AI 2 O 3 emis¬ 
sion coming from very close to the star and the super-solar abun¬ 
dances in the wind that are supposedly needed to account for the 
amorphous AI 2 O 3 emission, we consider in our model a gravita¬ 
tionally bound shell of amorphous AI 2 O 3 grains to be the main 
contributor to the emission at wavelengths around 12//m. We 
combine this gravitationally bound shell with an outflow to ob¬ 
tain a dust model that fits the wealth of information on the dusty 
wind of W Hya. The properties and spatial distribution of the 
dust in this model are discussed in the context of wind driving, 
as well as the time variability of the mass-loss on timescales of 
millennia. 


In Sect. 12 we detail our modelling strategy and the observa¬ 
tions of the dust envelope of W Hya. The modelling procedure 
and fits are described in Sect. [3] We discuss our results in Sect. 
|4]and provide a summary in Sect.|5] 


2. Observations of the dust envelope of W Hya 

In the following subsections, we present the dataset used in our 
study to constrain the present-day mass-loss rate. The observa¬ 
tions used to model the recent dust mass-loss history are pre¬ 
sented in Sect. 13.71 


2.1. Observed infrared spectra 


WHya was observed by the short-wavelength spectrometer 
(SWS, Ide Graau w et alJ 1996I) and the long-wavelength spec¬ 
trometer (LWS, Clegg etal. 1996I) o nboard the infrared space 
observatory (ISO, Kessler et alJll996h an d by the photodetector 
array camera and spectrometer (PACS, iPoglitsch et alj|2()10h 
and the spe c tral a nd photometric imaging receiver (SPIRE, 
|Griffin_etalJ_j20]^ onboard Herschel Space Observatory 
(iPilbratt et alJ l2010h . A spectrum covering the full range of 
the instruments was obtained with both PACS and SPIRE. The 
PACS observations were carried out in different epochs for the 
blue and red band. The full spectral range of SPIRE was covered 
in one observation run (see also Table[T]i. Two spectra were mea¬ 
sured with SWS and three with LWS. The two SWS spectra were 
taken almost one pulsation period apart (Tabled]). We obt ained 
the reduced 1997 spectrum from Llus ttanont e t al.l (l2004h and 
the 1996 spectrum from the ISloan et al.l (l2003h database. Two 
LWS observations were also taken roughly one period apart, 
with a third taken near maximum visual light. We obtained the 
highly processed data products of the three spectra from the ISO 
archive, in which the problem of near-infrared leakage is cor¬ 
rected. The spectrum obtained in Eebruary 1997 has an unreli¬ 
able baseline, especially at long wavelengths, and was not in¬ 
cluded in our analysis. The other two LWS spectra show very 
similar flux levels (displaying differences of about 10 %) and 
shapes. We have averaged these two spectra before fitting our 
models. All the observed spectra are shown in Eig. [1] the two 
SWS spectra differ slightly in flux level, revealing differences of 
up to 20%, but are very similar in shape. The PACS and SPIRE 
flux levels agree well with those of LWS when the different fields 
of view and the PACS far-infrared maps are taken into account 
(see Sect. [3. 

The dust components seen in the ISO spectrum of W Hya 
are better visualized if we subtract the expected stellar contin¬ 
uum from the observed spectrum. We used the stellar parame¬ 
ters and distance from Paper I, Tefj = 2500 K, L=5400Lo and 
78 pc, to estimate the stellar continuum. In Pig. l2l we compare 
the opacity curves of the d ust species considered bv ICamil (l2002h 
and iHeras & Honvl(l2005h and that of Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 to the 
two stellar-continuum-subtracted SWS spectra. A combination 
of amorphous silicates, amorphous AI 2 O 3 , and Mgo. 1 Peo. 9 O can 
qualitatively account for the emission seen between 9 pm and 
25 pm. Prom ~ 27 pm to ~ 32 pm, however, a feature is observed 
that cannot be reproduced by thes e dust species . We r efer to this 
feature as the broad 30 pm feature. Ipabi an et ^(l2001h identified 
a narrow feature at 32 pm that they assigned to spinel, but this is 
different from the broad hOpm feature we refer to here. 


2.2. Interferometric observations 

The envelope of W Hya was studied with the Masked Aperture- 
Plane Interference Telescope (MAPPIT) at the 3.9 m Anglo- 
A ustralian Telescope o ver the wavelength range 0.65 - 1.0jum 
by [Ireland et al.l (l2004h . The authors found an increase in mea¬ 
sured angular size from ~ 20 mas at 1 pm to ~ 35 mas at shorter 
wavelengths, which they proposed to be due the scattering of 
light b y dust grains. Later observations of W Hya bv iNorris et al.l 
( 12012 h using aperture-masked polarimetric interferometry with 
the NACO instrument on the VLT confirmed the existence of a 
shell of scattering dust particles around the star with a radius 
of about 37.5 mas, or ~2R*. The authors considered a popula¬ 
tion of large grains with uniform radii in their model, as their 
technique does not allow them to constrain the presence of small 
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Table 1. Overview of the observations of the dust shell of W Hya using different telescopes. The phase O of the pulsation period is 
given based on observations of visible light made available by the AAVSO“ 


Instrument 

Date 

0 

Ref. 

A [pm] 

FOV ["X"] 

Identifier 

ISO SWS 1 

14-02-1996 

0.12 

Sloan et al. 120031 

2.36-45.39 

14x20 to 17x40 

08902004 

ISO SWS 2 

07-01-1997 

0.17 

.lusttanont et al. 120041 

2.36-45.39 

14 X 20 to 17 x 40 

41800303 

ISO LWS 1 

07-02-1996 

0.1 

ISO data archive'’ 

43.05 - 195.5 

80 X 80 to 70 X 70 

08200208 

ISO LWS 2 

05-02-1997 

0.1 

ISO data archive’’ 

43.05 - 195.5 

80 X 80 to 70 X 70 

44700672 

ISO LWS 3 

02-08-1997 

0.6 

ISO data archive’’ 

43.05 - 195.5 

80 X 80 to 70 X 70 

62500572 

PACS Spec B2A 

09-07-2011 

0.8 

Khouri et al. 12014al 

55-73 

30x 30 

1342212604 

PACS Spec RIA 

09-07-2011 

0.8 

Khouri et al. 12014al 

102 - 146 

30x 30 

1342212604 

PACS Spec B2B 

14-01-2011 

0.4 

Khouri et al. 12014al 

70 - 105 

30x30 

1342223808 

PACS Spec RIB 

14-01-2011 

0.4 

Khouri et al. 12014al 

140-210 

30x30 

1342223808 

SPIRE Spec 

09-01-2010 

0.27 

Khouri et al. 12014al 

124-671 

37x37 to 18x 18 

1342189116 

PACS Image 70 

08-02-2011 

0.47 

Cox et al. 120121 

60-85 

1440 x 1440 

1342213848 

PACS Image 160 

08-02-2011 

0.47 

Cox et al. 120121 

130-210 

1440 x 1440 

1342213849 

SCUBA 




450 and 850 



VLT - MIDI 

04-2007 - 09-2009 

- 

Zhao-Geisler et al. 120111 

8-13 



VLT - NACO 

03-2009 and 06-2010 

0 .2" 

Norris et al. 120121 

1.04, 1.24 and 2.06 




Notes, a - Americ an Association of Variable Star Observers at http://www.aavso.org/, b - http://iso.esac.esa.int/iso/ida/, c - Visual phase given by 
iNorris et all(l2012h . 




Fig. 1. Compilation of all thermal emission observations of the 
dust in W Hya. The red and blue solid lines are the ISO SWS 
spectra from 1996 and 1997, respectively, the green solid line is 
the averaged ISO LWS spectrum, the purple and light blue solid 
lines are the PACS and SPIRE spectra, and the black squares and 
black circles depict the PACS and SCUBA photometric measure¬ 
ments. A black-body spectrum of 2500 K and 5400 Lq is over¬ 
plotted to guide the eye (black dashed line). 


Fig. 2. Continuum-subtracted ISO SWS spectra from 1996 and 
1997 are shown by the red and blue lines. The opacity curves 
of amorphous AI 2 O 3 , MgFeSi 04 , FeO, Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 
and Ca 2 Al 2 Si 07 are shown by the yellow, light brown, light 
blue, purple, and blue dashed lines, respectively. The opacity 
curves are scaled to compared them with the stellar-continuum- 
subtracted spectrum. 


grains, and found the scattering particles to be 0.316 + 0.004 jjm. 
Although the composition of the dust could not be determined, 
the authors suggested that they may be composed of iron-free sil¬ 
icates such as forsterite (Mg 2 Si 04 ) or enstatite (MgSiOs), while 
corundum (AI 2 O 3 ) is a possible candidate as well. 

IZhao-Geisler et al.l (1201 Ih monitored W Hya using MIDI at 
the VLT. The authors found that the silicate emission seen at 10 
yum in the ISO spectrum must come from a region with an in¬ 
ner radius larger than 500 mas, corresponding to 25 times the 
ste llar radius or 12 tim es the scattering shell radius as estimated 
bv iNorris et alJ(l2012h . This does not exclude the possibility that 
silicates condense close to th e star, as the iron-free silicates pro¬ 
posed bv iNorris et al.l(l2012h to be the scattering agent are trans¬ 
parent in the near infrared and can remain cold (~700K) and 
invisible even at two stellar radii. Such grains would not neces¬ 
sarily produce a signature in the thermal infrared (see Sect. 13.4b . 


2.3. Gas-phase wind model 

In Paper I, the present-day gas mass-loss rate of W Hya was de¬ 
termined to be (1.3+0.6)x10^^Mq year^' by fitting the strengths 
and shapes of the CO ground-vibrational state, pure-rotational 
transitions observed by Herschel. This mass-loss rate places con¬ 
straints on the amount of dust that can form in the wind assuming 
a solar composition. More stringent limits on the a mount of sil¬ 
icates t hat may form around W Hya were set by iKhouri et al.l 
(l2014bh . who modelled the SiO gas-phase emission lines ob¬ 
served by Herschel using the wind model obtained in Paper I. 
The authors found that about 60% of the silicon in the outflow 
of W Hya remains in the gas phase. 

3. Towards a comprehensive present-day 
mass-loss model 

Our modelling strategy consists of dealing with the different 
spatial and compositional components of the wind in an inde- 
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pendent way first. This approach can be followed as the wind 
of W Hya is optically thin at the wavelengths where dust emis¬ 
sion is important. We start by modelling the innermost observed 
dust component, namely the scattered light fractions and the 
amorphous AI 2 O 3 emission. Then, we model the silicate emis¬ 
sion envelope, whose inner radius is at ten times the distanc e 
from where the scattered light is seen (IZhao-Geisler et afl201 ih . 
Finally, we compare our best inner wind model with the PACS 
70/tm image, which provides constraints on both the dust dis- 
tribution close to the star and at scales up to a few arcminutes 
(ICox et alj|2012h . In our modelling efforts we aim to reproduce 
the following observations: 

- The s cattered light fractions as measured by iNorris et"^ 

dloTl . 

- The ISO dust spectrum. 

- The lower limit of 5 00 mas on the inner radius of the silicate 
emission shell set bv IZhao-Geisler et al.l (1201 Ih . 

- The limits set by t he gas-phase model of Paper I and 
iKhouri et alJ (l2014bl) on the amount of elements available 
for dust formation in the wind. 

- The broadness of the 70 ^m r adial brightness distribution 
seen by PACS (ICox et al.ll2012h . 


3.1. MCMax and model assumptions 


To solve the continuum radi ative transfer in th e wind of W Hya, 
we used the code MCMax (iMin et alj|2009h . The code calcu¬ 
lates the dust spectrum, images of the envelope at different wave¬ 
lengths, and the fraction of light scattered by the dust particles, 
all of which are used to constrain our model. The dust enve¬ 
lope was assumed to be spherically symmetric. To derive dust 
opacities from opti cal constants, we applied a distribution of hol¬ 
low spheres (DHS. lMin et al.ll2003h approximation to represent 
the particles shapes. We used the direction-dependent scatter¬ 
ing phase-function in our calculations, since particles with radii 
in the range considered scatter light non-isotropically. This was 
done by considering the full angle-dependent Mueller matrix for 
each scattering event. The temperature of the dust grains was set 
by assuming radiative equilibrium. For multiple dust species it 
may be computed assuming the different types of grains to be 
either in isolation or in thermal contact. In the case of thermal 
contact, the absorbed energy is distributed over the constituents 
of the grain and a single grain temperature is computed such 
that the emitted energy matches the absorbed energy. We spec¬ 
ify throughout the text whenever thermal contact is considered. 

To model the amorphous AI 2 O 3 emission, we used three 
sets of o ptical constants, those fo r porous and compact particles 
given by |Begemann_et^j2997l) and those of aerosil particles 
measured bv iKoike et al. (Il995h . The latter measurements were 


carried out for wavelengths between 0.5 and 500 u m. The op¬ 
tical constants obtained by iBegemann et al.l (Il997h were mea¬ 
sured between 7.7 and 500//m and lack data at near-infrared 
and optical wavelengths. As these short wavelength data are 
important for determining the tem perature of t he gra ins, we 
used the optical constants given bv iKoike et al.l (Il995h at d < 
7.7 pm. To model the silicate emission, we retrieved data for sil¬ 
icate dust species f rom the JENA dat a base of optical c o nstant s 
from the works of .lager et al. ( 1994h: iDorschner et al.l (Il995h : 
iMutschke et al.l ([19981) . and Jager et al.l ( 2003 ). The dust species 
that are considered are given in Table lA.ll For many of the sil¬ 
icate species the optical constants are available only down to 
about 6 pm. As most of the radiation that heats the grains is ab¬ 
sorbed at shorter wavelengths, we considered the near-infrared 



Fig. 3. Models considering full condensation of aluminum in a 
solar-composition outflow compared to the stellar-continuum- 
subtracted infrared excess (solid red). The aerosil AI 2 O 3 parti¬ 
cles, porous AI 2 O 3 particles and, compact AI 2 O 3 particles are 
represented by the dotted purple, short dashed blue and long 
dashed green lines. None of the considered species produces 
strong enough emission at 12 pm. 


optical constants of these dust species to be equal to that of 
Mgo. 4 Feo. 6 Si 03 , Mgo. 8 Feo. 2 Si 03 or MgSi 03 , selecting out of 
these three the species that best matches the iron content. 

We considered W Hya to be at 78 parsecs, to have a luminos¬ 
ity of 5400 Lq and the stellar spectrum to be that of a black-body 
with 2500 K, following our assumptions in Paper I. 

3.2. Previous dust model 

The model obtained in Paper I for the wind of WHya 
(shown in Fig. ad opted the dust components identified by 
lJusttanont et akl (l2004l) . In Paper I, the region of the ISO spec¬ 
trum between 8 and 30 pm of the ISO spectrum was fitted, but 
no attempt was ma de to fit the 13 um feature as its carri er is 
still uncertain (see iPosch et ^ Il999t IZeidler et al.l 1201 31 and 
references therein). A total dust mass-loss rate of 2.8 x 10“'° 
Mq year“* was obtaine d. The dust composition by mass is 58 % 
astronomical silicates (lJusttanont & Tielenall992h . 34 % amor¬ 
phous aluminum oxide (AI 2 O 3 ), and 8 % magnesium-iron oxide 
(MgFeO). Optical constants for the latter species were retrieved 
from the database of the University of Jena from the work of 
iHenning et al.l (Il995h . For amorphous aluminu m oxide, the au¬ 
thors h ave used data for porous particles from iBegemann et al.l 

dHH. 


3.3. Amorphous AI 2 O 3 emission 


As pointed out in the introduction, the amorphous AI 2 O 3 
mass fraction of the total dust needed to explain the mid- 
infrared e xcess of low mass-loss rate AGE stars is ofte n very 
high (e.g. Lorenz-Martins & Pomn^ 2000t ISneck et~^ 


iHeras & Honvl 2005 ; Karovicova et al. 20131) . The maximum 
amount of a given dust species that can exist in the stellar wind 
is determined by the availability of its least abundant element. 
In the case of amorphous AI 2 O 3 this limit is set by aluminum, 
whose mass abund ance relative to hydr ogen is 7.6 x 10“^ for a 
solar composition (lAsplund et al.lUOCM . If all aluminum atoms 
in the wind of WHya were used to form amorphous AI 2 O 3 
grains, the mass-loss rate of this dust species alone would be 
MAI 2 O 3 = 1.3 X 10“^ X 7.6x10-^0.53 = 1 . 9 xlO-"M 0 year-', 
where 0.53 is the fraction of the mass of amorphous AI 2 O 3 in 
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aluminum atoms and 1.3 x 10"^ Mq year^' the gas mass-loss 
rate. Assuming the specific weight of amorphous AI 2 O 3 to be 
4.0 g cm“^, we computed models with this amorphous AI 2 O 3 
mass loss, condensing at two stellar radii. At 12yum, the excess 
emission produced by these models is only 5 % of what is seen 
in the ISO spectrum (see Fig. [3]). These calculations depend on 
the adopted grain model, the optical constants used, and the as¬ 
sumed location of the onset of the wind. We investigated these 
parameters using the wind model derived in Paper I by modify¬ 
ing the opacities for amorphous AI 2 O 3 at short wavelengths and 
varying the inner radius of the dust envelope accordingly. We 
find that even an increase in the near-IR opacity of amorphous 
AI 2 O 3 by a factor of 10 increases the emission at 12//m by only 
about a factor of four. Increasing the near-IR opacity also leads 
to an important contribution of amorphous AI 2 O 3 emission at 
wavelengths shorter than 10/tm. We conclude that the emission 
attributed to amorphous AI 2 O 3 in W Hya cannot be accounted 
for by such grains in the present-day wind alone. 

3.4. Gravitationally bound dust shell 

We calculated models for a gravitationally bound dust shell 
(GBDS) located at about 2.0 R* to account for the amorphous 
AI 2 O 3 emission and the scattered light fractions. The GBDS is 
composed of grains located in the region between the radius of 
onset of condensation of amorphous AI 2 O 3 and that at which 
the radiation pressure on the dust grains overcomes the pull of 
gravity. Since we do not calculate the dust condensation nor the 
transfer of momentum to the dust grains, the inner and outer radii 
of this region are free parameters in our models. In the GBDS the 
grains must reach the optimal size for scattering, at which radi¬ 
ation forces can accelerate them. We therefore expect to find a 
distribution of sizes ranging from small grains, with radii of the 
order of 0.01 fim to larger grains with radii of about 0.3 //m. As 
an instructive exercise, we can approximate the size distribution 
by a mix of two particle sizes, small and large grains. These can 
be studied somewhat independently based on the observations 
available: large grains dominate the scattering of radiation, and 
small grains do not contribute to the scattered light and proba¬ 
bly dominate the infrared emission. We have considered small 
grains of 0.03 fim and large grains of 0.3 /tm. 

We considered the grains in the GBDS to be composed of 
either amorphous AI 2 O 3 or Mg 2 Si 04 . The amorphous AI 2 O 3 
emission can be fitted equally well with large and small grains 
of this dust species. Silicate emission from close to the star is 
not detected (Zhao-Geisler et al. 2011) but this does not neces¬ 
sarily mean that silicates are absent. Iron-free silicates (such as 
Mg 2 Si 04 ) can exist at close distances from the surface because 
they are translucent in the optical and near-infrared and therefore 
remain relatively cold. This does imply that we will not be able 
to constrain the mass in small Mg 2 Si 04 grains with our mod¬ 
els. A fit to the observed scattered light fractions requires large 
grains and in principle can be obtained with either type of grain. 
Such a fit sets an upper limit on the amount of large grains of the 
given species. 

We considered the amorphous AI 2 O 3 emission to originate 
from a stationary shell extending from 1.7 to 2.0 R*. The inner 
radius was set at that at which the small amorphous AI 2 O 3 grains 
reach temperatur es of 1400 K . The o uter radius was assumed to 
be that found by iNorris et~an ( 1201 2 h for the origin of the scat¬ 
tered light. In principle the grains in this region may be created 
and destroyed as material rises and falls as a consequence of stel¬ 
lar pulsations, but that does not affect our results. Our modelling 
shows that the amount of emission is only mildly sensitive to 


the exact value of the boundaries or to a radial density gradient 
within the shell, but mostly to the total mass of emitting amor¬ 
phous AI 2 O 3 . The scattered light fractions w ere fitted consider- 
ing a shell with an inner radius of 2.0 R* (iNorris et alJl201^ 
and thickness of 0.1 R*, that is, enclosing the shell from where 
amorphous AI 2 O 3 originates. In the context of our model, this 
shell of large grains may be understood as being the interface 
between the GBDS and the outflow. 

A population of only small amorphous AI 2 O 3 grains reaches 
temperatures of 1250K at 2.0 R*, while large amorphous AI 2 O 3 
grains reach 1500 K at the same distance. The higher temper¬ 
ature of large amorphous AI 2 O 3 grains has two causes. First, 
the mass-absorption coefficient between 0.2 and 2 fim is higher 
for the large grains than for small grains. Consequently, large 
grains absorb more radiation by mass than their smaller coun¬ 
terparts and reach higher temperatures. Second, the scattering of 
photons by large grains contributes to the diffuse radiation field 
and the density of short-wavelengths photons (A < 2 fim) sig¬ 
nificantly increases. The stronger radiation field at shorter wave¬ 
lengths leads to more energy absorbed per time per grain and to 
a higher dust temperature. This second effect also causes small 
amorphous AI 2 O 3 particles to have higher temperatures when 
placed in a region with 0.3//m grains. Assuming solar compo¬ 
sition and full aluminum condensation, a shell of small amor¬ 
phous AI 2 O 3 grains between 1.7 and 2.0 R* requires a gas mass 
of ~ 1.4 X 10“^ Mq to reproduce the amorphous AI 2 O 3 emis¬ 
sion seen in the ISO spectrum. Such a shell produces an infrared 
spectrum indistinguishable from the amorphous AI 2 O 3 models 
shown in Fig.|5]for a population with a range in grain radii. 

The mass in large grains needed to reproduce the scattering 
is 4.9 X 10“'° Mq for amorphous AI 2 O 3 and 3.6 x 10“'° Mq 
for Mg 2 Si 04 grains, assuming specific weights of 4 g cm“^ and 
3.2 g cm“^, respectively. By assuming a solar composition, full 
aluminum condensation, and 35% of silicon condensation (as 
to be consistent with the available sili con for dust condens ation 
derived from the gas-phase analysis; iKhouri et al.l 12014^ . we 
obtain gas masses of 3.5 x 10“° M© and 2.3 x 10“^ M© for large 
grains of amorphous AI 2 O 3 and Mg 2 Si 04 , respectively. 

If the fit to both the ISO infrared excess and the scattered 
light fractions are combined, the amorphous AI 2 O 3 emission 
and the scattered light fractions can be explained by a single 
population of amorphous AI 2 O 3 grains located in a shell with 
inner and outer radii of 1.7 and 2.1 R*, respectively, and with 
a grain size distribution of the type, n{a) oc where a is 

the grain radius. This expression corresponds to the standard 
distribution for grai n sizes found in the interstellar medium by 
iMathis et~^ (Il977l MRN distribution). The lower limit of the 
distribution was assumed to be 0.01 um a nd the maximum grain 
radius to be that derived bv INorris et al.L i.e. 0.316//m. By as¬ 
suming solar composition and full aluminum condensation, the 
total gas mass required for this population of amorphous AI 2 O 3 
is ~1.2 X 10“° M©. The differences in the amorphous AI 2 O 3 
mass needed to reproduce the observed emission only consider¬ 
ing small grains or a grain size distribution is due to the different 
grain temperatures obtained. The fit to the scattered light frac¬ 
tions using an MRN size distribution is very similar to the one 
obtained with the large-grains-only model, shown in Fig. |4] and 
is not given. Gradients in the grain size throughout the GBDS 
may be present and would allow similar quality fits as adopting a 
size distribution. Infrared spectra for an MRN distribution of par¬ 
ticles adopting different optical constants for amorphous AI 2 O 3 
are shown in Fig. 5. The similar amorphous AI 2 O 3 emission ob¬ 
tained shows that the applied optical constants do not impact our 
findings significantly either. 
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Fig. 4. Scattering fractions for the models with large particles of 
aerosil AI 2 O 3 (red long dashed line) and Mg 2 Si 04 (blue short 
dashed line) in the GBPS. The black squares are the observa¬ 
tions bv iNorris et al.l (l2012h . 
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If the scattering is due to large Mg 2 Si 04 grains, a population 
of small amorphous AI 2 O 3 grains is still required to account for 
the infrared excess emission. When a population of large sili¬ 
cate grains is placed in the shell together with small amorphous 
AI 2 O 3 particles, the temperatures are about TOOK and 1400K, 
respectively. When these two species are placed in thermal con¬ 
tact, the resulting temperature will depend on the mass fraction 
of each species in a given grain. Assuming full aluminum and 
35% silicon condensation implies that 8 % of the mass is in amor¬ 
phous AI 2 O 3 and 92% in Mg 2 Si 04 . In this case, the equilibrium 
temperature is about 800 K and the silicate features would re¬ 
main undetected, in agreement with the limit set for the inner 
radius of amorphous silicate emission of ~40AU. Therefore, 
a scenario in which small amorphous AI 2 O 3 grains are formed 
close to the star and in which rapid condensation of Mg 2 Si 04 
onto the amorphous AI 2 O 3 seeds occurs as well would be con¬ 
sistent with all observations. In this scenario, the condensation 
of Mg 2 Si 04 could be halted by the onset of the wind. 



Fig. 5. Spectra of a population of amorphous AI 2 O 3 , MgSi 03 
or Mg 2 Si 04 grain s compared to the observed ISO spectrum 
(ISloan et al.l l2003h . We considered an MRN grain radii dis¬ 
tribution for all models shown here. The different lines are 
for grains consisting of aerosil AI 2 O 3 particles (green dashed), 
AI 2 O 3 porous particles (light blue long dashed-dotted), AI 2 O 3 
compact particles (brown short dashed-dotted), MgSi 03 (purple 
dotted) and Mg 2 Si 04 (short dashed blue). 


3.5. Dust properties in the outflow 

The contribution of amorphous AI 2 O 3 grains in a GBDS to the 
infrared spectrum is somewhat independent of the size distribu¬ 
tion and the exact location of the shell as the grains would be 
sufficiently warm for the shape of their emission feature to be 
independent of temperature for the optical constant adopted by 
us. Although the opacities o f crystalline AI 2 O 3 s how significant 
variations with temperature (IZeidler et al.ll2013h . to our knowl¬ 
edge, no studies have been performed that probe the effect of 
temperature variations on the amorphous AI 2 O 3 features. In the 
context of our model, only the strength of the emission is sensi¬ 
tive to temperature, and therefore the total grain mass needed to 
reproduce the infrared excess is a function of grain temperature. 
If we assume that the dust shell is optically thin, we may sub¬ 
tract the amorphous AI 2 O 3 contribution from the ISO spectrum 
and study the residual spectrum (see Fig. |6ll. The assumed size 
distribution of amorphous AI 2 O 3 grains does not affect the over¬ 
all shape of the residual spectrum. Uncertainties in the residual 
spectrum arise from the choice of optical constants for amor¬ 
phous AI 2 O 3 , the adopted grain model, and variability of the star, 
but these probably have a modest effect. 

The residual spectrum clearly shows the silicate peak at 9.7 
pm, the unidentified 13 pm feature, a peak at about 20 pm that 
is probably due to 18 pm silicate and/or magnesium-iron ox¬ 
ide emission, and a broad feature around 30 pm. For simplicity 
and because of the similarity in the different Al 203 -subtracted 
resi dual spectra , we on ly used the ISO spectrum retrieved from 
the ISloan et ^ (l2003h database subtracted by the spectrum of 
aerosil amorphous AI 2 O 3 grains (Fig.| 6 j red solid line). 


3.5.1. Fit to the residual spectrum 

The dust species considered in fitting the silicate emission are 
given in Table lA.ll All the species listed have somewhat sim¬ 
ilar opacity profiles that are characteristic of silicate glasses. 
Diagnostic parameters such as the peak strength ratio, the central 
wavelength, and the full-width at half ma ximum (FWHM) of the 
peaks can help constra in the composition (iDorschner et al.ll995l: 
iMutschke et al.l Il998h and are listed for DHS particles with 
0.3/im radius. The position of the two silicate peaks measured 
from the Al 203 -subtracted residual spectra are 9.94 + 0.5 pm and 
19.7 + 0.5 pm. It is important to note that the peak at 19.7 pm 
might be produced by a combination of silicate and MgxFei_xO 
emission, so the measured position should be considered with 
care. 

We recall that the VLT 10 pm observations provide a lower 
limit to the radial distance of the silicate-emitting shell of about 
40 AU. Most of the dust species listed in Table lA.ll have a 
low iron content and consequently a low near-infrared opacity. 
These species would have too low temperatures at such large dis¬ 
tances and would not produce detectable infrared flux relative to 
the bright stellar continuum and emission from the other dust 
species. As the ratio between the two silicate emission peaks 
indicates that the grains must be warm (~500K), the silicates 
with low near-infrared opacities, that is, with low iron content, 
were placed in thermal contact with metallic iron to increase 
their temperatures. Such metallic iron inclusi ons are indeed ex¬ 
pecte d to form in circumstellar silicates grains (iGail & Sedlmavil 
Il999h . We considered models with an iron content per unit mass 
ranging from 0 to 80%. Accounting for even higher amounts 
of metallic iron causes the silicate features to become too weak 
with respect to the continuum emission that the metallic iron 
produces. Silicate species with significant iron content, such as 
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Fig. 6. Al 203 -subtracted residual infrared spectrum of W Hya. 
Upper panel: Residuals adopting the optical constants for aerosil 
particles fromiKoike et al.l (Il995h and t he ISO spectra from the 
ISloan et ^ (l2003ll database (in red) and Llusttanont et all (l2004ll 
(in blue); lower panel: The red curve is the same as in the upper 
panel. The yellow and blue curve s are for residua l s assu ming 
compact and porous particles from iBegemann et al.l (Il997h . 

MgFeSi 04 , have high enough temperatures even when placed at 
40 AU and require modest or no thermal contact with metallic 
iron at all. 

Models were calculated for silicate envelopes ranging from 
40 AU to 1800 AU (i.e. the outer limit given by the field of 
view of ISO SWS) for all species listed in Table I A. ll separately. 
When in thermal contact with substantial amounts of metallic 
iron, the silicate grains can be heated up to 400 - 450 K at 40 
AU and some of the species considered produce peak ratios that 
reflect the observed ratio. As the temperatures that are reached 
are still fairly low, species that have relatively high values for 
the are preferred. Furthermore, metallic iron pro¬ 

duces emission that fills the trough between the silicate peaks 
and matches the observed AUOs-subtracted residual spectrum. 
Figure I?] shows the effect of varying the amount of iron in ther¬ 
mal contact with the silicate grains. A higher iron content im¬ 
plies more emission in the silicate trough and blue-wards of 
the 9.7 pm peak. For a given silicate species, the ratio between 
the red and the blue silicate peaks decreases with an increasing 
amount of metallic iron inclusions. The dependence is stronger 
for smaller amounts of iron, between 0 and 50% per mass. 

When using a continuous distribution of ellipsoids (CDE, 
iBohren & HuffmanI 119981: iMin et all l2003h as the grain model 
for the metallic iron particles, the amount of this species 
needed to heat the silicates to a given temperature is greatly re- 
duced if compared t o models with DHS metallic iron particles. 
iKemner et al.l (l2002h also found a large difference between the 


opacity of CDE metallic iron grains and that obtained using Mie 
theory. We find the DHS metallic iron model to have an interme¬ 
diate opacity between these two extremes and discuss the effects 
of using CDE metallic iron grains in Sect. 14.3.11 



Fig. 7. Models for an envelope composed of Mg 2 AlSi 207.5 and 
metallic iron ranging from 40 AU to 1800 AU. The different lines 
are for 30% (dotted purple), 50% (short dashed blue) and 70% 
(long dashed green) metallic iron content. The AI 2 O 3 -subtracted 
residual ISO spectrum is shown by the full red line. 

Another feature in the Al 203 -subtracted residual spectrum 
worth noting is the broad component at 30 pm. The mass-loss 
rate of W Hya varies on larger scales, therefore this emission 
bump might be thought to be a signature of an abrupt change in 
mass-loss rate. However, the bump seems too sharp to be pro¬ 
duced in such a way. Moreover, this feature is seen in many 
other ISO spectra of oxygen-rich AGB stars wi th low mass- 
loss ra te, suggesting a more fundamental cause. iGervais et al.l 
(Il987h studied the infrared properties of silicate glasses and con¬ 
cluded that calcium-bearing silic ates have a characteris tic band 
around 30 pm. From the work of iMutschke et al.l (Il998h . we re¬ 
trieved the optical data for two Ca-silicates, which do indeed 
show a feature at 30 pm. Ca-bearing silicates are expected to 
form from gas-solid chemical rea ctions between AI2O3 grain s 
and calcium and silicon in the gas (iGrossman & Larimeilll974ll . 
Another candidate for explaining the feature is pure calcium- 
oxide grains. However, this species produces a rather localized 
peak, whereas calcium-bearing silicates also provide continuum 
opacity beyond 30 pm, which better matches what is seen in the 
spectrum of W Hya. 

As shown in Fig. [ 8 j the calcium-aluminum silicate fits the 
Al 203 -subtracted residual spectral shape nicely when in ther¬ 
mal contact with metallic iron. Unfortunately, we were unable 
to find other optical constants of calcium-bearing silicates in 
the literature, which potentially could provide a better fit to the 
9.7 pm silicate peak. As calcium is indeed expected to condense 
just after AI 2 O 3 , the existence of this dust species in the wind 
of W Hya can be expected. The model shown in Fig. [ 8 ] has a 
dust mass-loss rate of 3 x 10“'° Moyear^', consisting of 60% 
Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 and 40% metallic iron, per mass. Further 
constraints on these quantities are obtained in Sect. l3.6l bv fitting 
the broadness of the PACS IQ pm brightness profile. 

Even though the fit using Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 is offset in the 
9.7 pm region, this single dust species is able to qualitatively 
reproduce the Al 203 -subtracted residual spectrum very well. 
Therefore we conclude that Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 is the silicate 
from our set of optical constants that best represents the com¬ 
position of the silicate dust in the envelope of W Hya. A better 
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Fig. 8 . Fit (dashed blue) to the amorphous Al 203 -shell- 
subtracted thermal spectrum (in solid red), adopting 40 AU as 
the inner boundary of the silicate emitting envelope and a mix¬ 
ture of 40% metallic iron and 60% Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 particles 
that are in thermal contact. 


fit to the individual features might require a more complex dust 
envelope structure, which could be linked to the condensation of 
different elements at different distances from the star, but such 
a more complex dust model is beyond the scope of this paper. 
Having this amount of Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 condense from the 
outflowing material, however, requires a super-solar abundance 
of calcium and aluminum. We discuss this in Sect. I4.3.2I 

3.6. Fitting the peak vaiue and broadness of the PACS 70pm 
image 

When the model consisting of an amorphous AI 2 O 3 GBDS and 
Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 grains in the wind that are in thermal con¬ 
tact with metallic iron (see Fi g. is compar ed to the PACS 
70 pm radial brightness profile dCox et ani2012h . the peak value 
of this profile is not reproduced. The model peak flux is 30% 
lower than the observed value and essentially constrains the 
present-day dust mass-loss. The broadness of the intensity pro¬ 
file in the first few arcseconds constrains the distribution of the 
bulk of the 70 pm emission. 

When matching the present-day dust mass loss to the peak 
strength, the fit to the ISO spectrum becomes poorer. To recover 
the fit, we shifted the inner radius of the silicate shell from 40 
AU to 50 AU and increased the iron content from 40% to 50% 
of the dust mass. The width of the PACS intensity profiles is 
matched if we decrease the dust mass-loss rate beyond 500 AU 
by a factor of two, to 2 x IO^'^Mq year^*, preserving the dust 
composition from the inner part. The best model thus consists of 
a wind composed of 50% Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 and 50% metallic 
iron in thermal contact and mass-loss rates of 4x10“'° M 0 year“' 
between 50 AU and 500 AU and of 2 x 10“'° Mq year“' from 500 
AU to 1800 AU. The fit to the ISO SWS spectrum of this best 
model for the outflow combined with the GBDS of amorphous 
AI 2 O 3 is shown in Fig. [10] Comparing our dust model and the 
CO model (Paper I) yields a gas-to-dust ratio in the inner wind 
of W Hya of 325. 

3.6.1. CO envelope 

The mass-loss rate discontinuity at 500 AU is close enough to 
the star that it might affect the CO envelope. The model for the 
CO rotational emission presented in Paper I does not take this 
variable mass loss into account and therefore has to be revised. In 


Fig. 9. Fit to the radial intensity profiles observed by PACS 
at 70 pm. The blue and green solid lines represent the ob¬ 
served brigh tness distr i bution averaged over a narrow range 
of direction (ICox et al.ll2012h and over all directions, respec¬ 
tively. The long dashed red line represents the model with M - 
3 X 10“'° Mq year“', inner and outer radius of 40 AU and 1800 
AU, respectively, and with grains composed of 40% metallic 
iron and 60% Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 . The short dashed black line 
shows the model with M — 4x 10“'° Mq year“', inner and outer 
radius of 50 AU and 1800 AU, respectively, and with grains 
composed of 50% metallic iron and 50% Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 . 
The dotted purple line represents the model for the recent mass- 
loss history obtained in Sect. 13.71 with parameters given in Fig. 
M and with grains composed of 50% metallic iron and 50% 
Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 . The radial distance in AU is given along 
the upper x-axis for an assumed distance of 78 parsecs. 



Fig. 10. Overall fit of the best model (short dashed dark blue line) 
to the ISO SWS spectrum (solid red line). The model consists of 
the GBDS composed of amorphous AI 2 O 3 grains plus an outflow 
(50% Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 and 50% metallic iron) in which the 
mass-loss rate is 4x 10 “'°Mq year“' from 50 AU to 500 AU and 
2 X 10“'°Mo year“' from 500 AU to 1800 AU. The stellar input 
black-body spectrum is shown by a dashed light blue line. For 
comparison the best model from Paper I is shown by the long 
dashed green line. 


the CO study a fit to both the high- and low-excitation rotational 
lines was achieved by decreasing the photodissociation radius of 
CO by a factor of 2.5, placing it at about 500 AU from W Hya. 

In Fig. [TT] we compare the new model for the CO emission 
that adopts the new dust model and the jump in the gas-mass- 
loss rate from 1.3 x 10“^ to 7 x 10“®Moyear“' at 500 AU. The 
new CO model reproduces the CO rotational lines as well as a 
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Fig. 12. Schematic dust mass-loss history of W Hya. The mass- 
loss rate shown by the dashed line could not be constrained by 
our model. The timescales given along the upper x-axis is calcu¬ 
lated assuming a constant expansion velocity of 7.5 km s“'. 


Fig. 11. We compare the observed CO rotational emission lines 
(red solid line) to the model with mass-loss properties as de¬ 
scribed in the caption of Fig. fTOland a dissociation r adius set to 
the canonical value obtained bv iMamon et al.l (Il988h (shown by 
the blue dotted line). The model obtained in Paper I is shown for 
comparison (black dashed line). 


model with CO dissociation occurring at 500 AU and thus is in 
concordance with such a M discontinuity at this distance. 

3.7. Recent (dust) mass-loss history of WHya 

The mass-loss history over the past 10^-10^ years is unveiled 
by the infrared images obtained with PACS and SPIRE and th e 
infrared astronomical satellite (IRAS, iNeugebauer et al.lll984h . 
The image taken by PACS at IQfim has the best spatial res¬ 
olution and, therefore, provides the best constraints. Emission 
from the extended dust envelope that seems to be produced by 
a highe r mass-loss rate phase is also seen in the infrared ISO 
spectra (lHawkinslll990h . This emission is only produced by cold 
(T < 70 K) dust and affects the flux measured at the relatively 
long wavelengths (A > 45 fim) of the LWS spectral range; it has 
no impact on the dust emission in the inner wind. 

3.7.1. Infrared images 

The 70 pm PACS brig htness distributio n is rather unusual for an 
oxygen-rich AGB star (ICox et al.ll2012l) . These images show ex¬ 
cess dust emission produced between 30" and 100" that cannot 
be understood in terms of a constant mass-loss wind. In Fig.0we 
show the angle-averaged radial profile of the observed flux den¬ 
sity at 70 pm. The PACS image shows the envelope to have an 
ellipsoidal shape when projected on the sky. This causes the ra¬ 
dial profile averaged over all directions to have a shallower slope 
from around 75" on than radial profiles computed for a narrow 
range of directions. The difference in continuum level of about 
20% between PACS and ISO at 70 pm is an effect of the signifi¬ 
cantly smaller field of view of the PACS spectrometer compared 
to that of ISO LWS. 

On a larger scale, 100 pm maps obtained with IRAS show 
a very extended dust shell, 30' in diameter, w hich requires a 
high mass-loss rate to explain it (lHawkinslll990h . Together, the 


PACS and IRAS infrared maps show that the dust mass-loss rate 
of W Hya has not been constant in the past 10"^ to 10^ years. 

3.7.2. Model for the recent mass-loss history 

To model this extended emission, we used the inner wind dust 
composition and applied variations in the dust mass-loss rate. 
Since the PACS image shows an elliptical dust envelope and our 
model assumes a spherically symmetric wind, we have measured 
the radial brightness profile in a narrow range of directions (from 
130° to 170° measured from north to east) and used that to con¬ 
strain our model. Measuring the radial profile in a different di¬ 
rection produces a similar curve with the peak seen at 75" in Fig. 
0 shifted to larger distances from the star. The chosen direction 
corresponds to the direction in which the radius of the observed 
circumstellar shell is the smallest. The asymmetry seen in the 
PACS image could, in principle, be due to a direction dependent 
expansion velocity, but that is not considered in our modelling. 
The largest radius measured is larger by about 25% than the one 
modelled by us. A difference of 25% in radius translates into a 
difference in expansion velocity of also 25%. 

The dotted purple line in Fig. |9] shows the fit to the pro¬ 
file composed of a combination of the constant mass-loss rate 
epochs shown in Fig. [12] The mass-loss history and the con¬ 
straints motivating this model are given in Table |2 The values 
given assume the present-day expansion velocity of 7.5kms“', 
measured from the gas lines. If the gas-to-dust ratio remained 
constant (at 325), the dust mass-loss rates of phases II, III, and 
IV would correspond to 6.5 x 10^*, 3.9 x 10'^, and 2.0 x 10~® 
Mq year^' of gas. The exact value of the derived dust mass-loss 
rates depend on the model for the present-day dust mass-loss. If 
the composition of the dust producing emission at 5000 AU 
is not significantly different from that of the inner silicate en¬ 
velope (r > 50 AU), the derived ratio between the high and 
low mass-loss rate phases should be accurate even if our inner- 
wind model were not. If the dust composition and/or expansion 
velocity changes significantly between these two components, 
our results probably underestimate the mass-loss rate for the 
phase of stronger mass loss. That is because, first, the opacity 
of Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 at 70//m is higher than that of other sili¬ 
cates, so more mass of other silicate species at the same temper¬ 
ature would be required to produce the same emission at these 
wavelengths. And, second, the expansion veloc ity tends to in- 
crease with increasing mass-loss rate (see, e.g.. lOlofsson et aP 
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Table 2. Mass-loss history of W Hya. 


Phase 

Mdu.st 

[Mgyear^'] 

Rin 

[AU] 

Constraint 

I 

4.0 X 

50 

ISO spectrum 

II 

2.0 X 10-'“ 

500 

70 gm profile peak - 1 - CO model 

III 

1.2 X 10-“ 

4000 

70 gm profile from 20 to 60 " 

IV 

6.0 X 10-“ 

5400 

70 gm profile beyond 60 " 

V 

1.2 X 10-“ 

5900 

70 gm profile beyond 60 " 


I 2 OO 2 I) . which means that the present value of 7.5 km s ' would 
probably be lower than the expansion velocity during a phase 
with higher mass-loss rate. 

4. Discussion 

We have modelled the dust envelope of W Hya using two com¬ 
ponents: a GBDS from which the AI 2 O 3 emission originates and 
an outflow that contains silicates that are in thermal contact with 
metallic iron. Overall, the best-fit model (see Fig. 10) provides 
a good fit to the infrared emission, although the position of the 
silicate peak at 9.7 jim is not well reproduced. Our model does 
not require a super-solar aluminum abundance that is required if 
the AI 2 O 3 grains that cause the emission at 12/zm are part of the 
outflow. The silicate emission originates from beyond 50 AU and 
best fits the spectrum if silicate grains that contain Al and Ca are 
used. This again introduces a super-solar abundance problem but 
we argue that this may be a spurious problem. We now discuss 
the empirically derived envelope structure in the context of wind 
dynamics, the gas- and solid-phase silicon budget and the ele¬ 
mental abundances of Al and Ca, and the expected crystallinity 
of the dust grains. 

4.1. Scattering agents and wind driving in oxygen-rich AGB 
stars with a iow mass-ioss rate 

The spectrum of W Hya is well-explained by a GBDS of amor¬ 
phous AI 2 O 3 that produces most of the thermal emission seen 
at 12 yum. For the fit to the scattered light fractions our mod¬ 
els can not distinguish between the two species that have been 
suggested, namely silicates or AI 2 O 3 . Although AI 2 O 3 is able 
to explain the scattered light and thermal emission simultane¬ 
ously, wind-driving models suggest that a population of pure 
amorphous AI 2 O 3 grains provides insufficient opacity to initi¬ 
ate a wind. Silicates, however, are much more effective in ini- 
tiating the outflow bec ause of the higher abundance of silicon 
(iBladh & Hofneijl2012l) . To illustrate this: condensation of 35% 
of the available silicon in a solar composition gas into large 
grains of a typical silicate glass increases the scattering opac¬ 
ity by an order of magnitude relative to dust that is composed 
of fully condensed large amorphous AI 2 O 3 grains alone. If large 
amorphous AI 2 O 3 grains are lacking in the GBDS the 12 gm 
emission is produced mainly or exclusively by small amorphous 
AI 2 O 3 grains. 

The gas mass that should be present in the GBDS to ac¬ 
count for the amorphous AI 2 O 3 emission is ~ 10“^ M©. Given 
the mass-loss rate of WHya of 1.3 x 10“^ Mq year^*, this im¬ 
plies that only a small fraction (a few percent) of the GBDS is 
expelled per pulsation cycle of about a year. The gas mass in 
the GBDS that is required to produce the scattered light - as¬ 
suming it is due to 35% silicon condensing in large Mg 2 Si 04 
grains - is ~2.3 x 10^^ Mq. This mass is more than an order 
of magnitude lower than the corresponding mass for amorphous 


AI 2 O 3 , and is of about the same magnitude as the total mass that 
is lost per year. If the silicates are responsible for the scatter¬ 
ing and condense on the amorphous AI 2 O 3 seeds, condensation 
of this material should occur over a small radial interval of the 
GBDS. This is because the composite grains must be cold not to 
produce silicate emission and, therefore, would also not produce 
any amorphous AI 2 O 3 emission. Because of this and as a result 
of the low abundance of aluminum, a relatively larger fraction 
of the GBDS volume in which warm amorphous AI 2 O 3 grains 
exist is needed to reproduce the observed emission. Models of 
outflows driven by large Mg 2 Si 04 particles indeed show very 
rapi d growth of th ese grains in the onset region of the flow (see, 
e.g.. lHofn^l2008h . 

Alternatively, the grains responsible for the scattering could 
be aggregates of different oxides including species other than 
silicates and AI 2 O 3 - such as calcium or magnesium oxides - 
but we have not explored this possibility. 

4.2. Structure of the dust enveiope of WHya 

A schematic view of the model for the envelope of WHya 
is shown in Fig. [13] Most of the amorphous AI 2 O 3 emission 
can be accounted for by grains in the GBDS. The thickness of 
this GBDS is difficult to assess as its inner radius depends on 
the poorly constrained near-infrared optical constants of amor¬ 
phous AI 2 O 3 . At about 2R* grains of 0.3//m are present to 
account for the scattered light fractions observed. In our mod¬ 
els, these grains can be either amorphous AI 2 O 3 or iron-free 
silicates. If these scatterers are silicates, it is likely that the 
amo rphous AI 2 O 3 grains serve as seeds for silicate condensa¬ 
tion dKozasa & Sogawalll997bh . If silicates exist so close to the 
star, they have to be translucent in the near-infrared and remain 
cold such as not to produce significant emission. Combined with 
the broadness of the brightness peak in the PACS 70/rm image, 
this implies that the bulk of silicate emission at 9.7 gm has to 
come from an envelope with an inner radius of 50 AU. The gas 
temperature at this distance is ~350K dKhouri et alJl2014H) . 

If the silicates are the scatterers and are present near the 
star without producing significant emission, there would have 
to be a physical process that increases their temperature beyond 
40 AU. Such a process could be a change in their lattice struc¬ 
ture, which increases the absorption in the near-infrared, or the 
establishing of contact with a more opaque dust species, for in¬ 
stance small iron grains - creating a silicate with iron inclusions. 
Alternatively, if the scattering is caused by AI 2 O 3 , or grains com¬ 
posed of a combination of other oxides, the inner radius of the 
silicate emission envelope could be set by silicate condensation. 
It is not clear, however, why any of the processes mentioned 
would occur at such large distances. This interpretation is based 
on the assumption of a steady-state wind model for the inner 
regions, r < 1000 AU. The observed radial dependency of the 
dust-envelope parameters of W Hya could potentially also be the 
result of variable conditions at the wind acceleration region, but 
this explanation seems less likely (see Sect. lO l. 

4.3. Bringing the dust and gas modeis together: eiementai 
abundances in the inner wind. 

The abundances of the elemental species that are required to re¬ 
produce the thermal emission from the wind of W Hya are 4.8 
times the solar-composition value for calcium, 0.06 for magne¬ 
sium, 3.5 for aluminum, 0.24 for silicon, 0.07 for oxygen, and 
1.15 for iron. If the opacity of the metallic iron grains is com- 


10 
















Khouri et al.: Dusty wind of W Hya. 


, STAR 
/ , GBDS 

^ / SCATTERING SHELL 

8r 


' DUST-FREE 


NO SILICATE 

SILICATE 

EMISSION / 

EMISSION 


Fig. 13. Schematic view of the envelope of W Hya. The radii of 
the different structures are given to scale. The radius of the star 
and the scattering shell are 1.8 AU and 3.6 AU. The inner radius 
of the GBDS is not constrained by our model and is assumed to 
be 3.1 AU (1.7R*) for the small amorphous AI 2 O 3 grains. The 
inner radius of the silicate envelope show n is the lower limit of 
40 AU given by IZhao-Geisler et all (1201 Ih . 


puted using a CDE distribution of sizes, a factor of three less 
iron is needed. In what follows, we discuss the abundances of 
the dust species in more detail and offer potential solutions to 
the apparent super-solar abundances of calcium and aluminum. 

4.3.1. Silicon abundance 

I Khouri et al.l (l2014bh modelled the gas-phase SiO emission and 
found that between 15 and 50% of the silicon atoms are expected 
to be depleted from the gas-phase. When considered together 
with our dust model this accounts for all the silicon expected for 
a solar composition wind. 

However, since the calcium-bearing silicate that we adopted 
(Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 ) does not provide a perfect fit to the 9.7 yum 
silicate peak, the constraints on the solid-phase silicon abun¬ 
dance are not stringent. To assess this problem, we have com¬ 
pared models of other silicates that provide a better fit to the 
9.7 peak silicate feature to the observed spectrum. These dust 
species provide a poorer fit to the overall spectrum, however, as 
they reproduce neither the ISyum silicate peak nor the 30 yum 
feature. We have chosen to fit the flux ratio between the 9.7 
silicate peak and the trough at 11.9 and 14.5 //m. The trough 
flux is sensitive to the amount of iron in thermal contact with the 
silicates, which needs to be 50% per mass, typically, in our mod¬ 
els. The best fits are always achieved with in between 30% and 
40% of the silicon available from the gas phase. Therefore, the 
amount of silicon needed in the dust ag rees with the depletion 
fraction reported bv iKhouri et alJ (12014bh independent of the ex¬ 
act composition of the silicate. This is in virtue of the strength 
of the 9.7 peak correla ting well with the Si 02 co ntent of the 
given silicate species ('e.g. lKoike & HasegawalfT98^ . 

4.3.2. Aluminum and calcium abundances 

In our model, aluminum is condensed in amorphous AI 2 O 3 
grains in the GBDS and also in silicates present in the wind. 
We do not require amorphous AI 2 O 3 to be present in the wind 
but a significant amount in terms of aluminum abundance could 
be driven out without producing a strong spectral signature. We 
considered only the outflowing aluminum-bearing grains in the 
aluminum budget, as the grains that form in the GBDS form 
from a gravitationally bound gas mass that is not independently 
constrained by our models. Nonetheless, our dust model requires 


3.5 times more aluminum than what is available from the gas 
phase. Unfortunately, we do not have aluminum-free calcium¬ 
bearing silicates in our dataset of optical constants, so we cannot 
directly probe the impact of the alu minum content on th e optical 
properties of Ca 7 Mgn sAUSii sO?. iGervais et al.l (Il987h studied 
the influence of aluminum content on silicate emission features 
and found it to have a minor impact on the strongest ressonances. 
This suggests that the over-abundance of aluminum found by 
us may be spurious and that an equally good fit to the AI 2 O 3 - 
subtracted residual spectrum could be achieved with a silicate 
with a lower aluminum content. 

The required super-solar abundance of calcium is some¬ 
what more intriguing, as this element seems to have a signifi¬ 
cant impact on the strengths and shapes of the silicate features. 
Particularly the 30 //m feature is linked to the presence of this 
element. Like in the case of aluminum, however, we do not have 
a library of optical constants for silicate species with different 
amounts of calcium content, which prevents us from assessing 
its relation to the 30 yum feature. 

We identify three potential solutions to the problem of the 
super-solar abundance of aluminum and calcium in the dust 
found by us. The first option is that silicates that contain smaller 
amounts of cal cium and aluminum are responsible for the ther¬ 
mal emission. ISneck et^ (1201 ih presented measurements of 
silicates with a range of compositions and calcium and alu¬ 
minum content. One of the species studied by them, named 
basalt in their work, has calcium-to-silicon and aluminum-to- 
silicon ratios very similar to those expected for 20 % silicon con¬ 
densation and full calcium and aluminum condensation. The ab¬ 
sorbance spectrum of this species shown in Fig. 4 of their paper 
peaks at 10 jum, which is closer to the observed peak position 
for W Hya than that of Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 . The 18 yum peak of 
this species is also in accordance with the observed spectrum of 
WHya. Unfortunately for our purposes, the authors focused on 
analysing the 10 yum and 18 yum region and we have no infor¬ 
mation on the spectra of this species at 30 yum. Measurements of 
optical data for silicates with these characteristics are needed to 
determine whether the spectrum of W Hya can be reproduced by 
such a species. 

A second possibility to explain the apparent super-solar 
abundance of Al and Ca is that some of the features fitted by 
Ca 2 Mgo. 5 Al 2 Si 1 . 5 O 7 , such as those at 20yum and 30yum, are par¬ 
tially produced in the GBDS, similarly to AI 2 O 3 . In that case 
the abundance problem would not exist. This can be the case if 
other oxides are able to condense in the GBDS. Calcium is in¬ 
deed expected to be one of the fi rst elements to react whh AI 7 O 3 
from condensation calculations (iGrossman & Larimeiill974h . If 
calcium reacts with and is incorporated into the AI 2 O 3 grains be¬ 
fore the onset of the wind, other features in the spectrum could 
also be produced in the extended atmosphere. Furthermore, the 
addition of calcium to the dust grains in the GBDS could help 
driving the wind. However, the increase in the dust-to-gas ratio 
due to calcium condensation would only be of about a factor of 
two. 

Finally, calcium and aluminum could be overabundant in the 
wind if these atoms are included in the dust grains before the 
wind is fully driven and if only a fraction of the gas is driven 
out by these grains. Although this might be the case, this sce¬ 
nario seems unlikely as the dust grains should be well coupled 
to the gas at the densities expected in the region where the flow 
is initiated. 
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4.4. Lattice structure of the grains in the GBDS 

The grains in the GBDS must form from high-temperature gas 
and are heated by the stellar radiation field to temperatures that 
may differ from that of the gas. Being subjected to high enough 
temperatures with respect to the glass temperature of the given 
dust species will cause the dust grains to become crystalline. 
Whether the lattice structures of the dust grains are amorphous 
or crystalline is directly connected to the features produced by 
them in the infrared spectrum. In particular, our model requires 
that amorphous AI 2 O 3 grains exist close to the star (~ 2R*) to 
reproduce the ISO spectrum. As we obtain the grain temperature 
as a function of radius, we can speculate what the lattice struc¬ 
ture of AI 2 O 3 and silicate grains in the GBDS would be and test 
if that is consistent with our assumptions. 

4.4.1. AI 2 O 3 grains 

The physical conditions in AGB outfl ows favour the conde nsa- 
tion of amorphous AI 2 O 3 grains (e.g. iDeir Adi et alllMil . In 
our calculations the temperature of the grains in the GBDS is 
found to be between 1500 K and 1600 K for large grains and 
betwe en 1200 K and 1300 K for small grains. iBegemann et alJ 
(Il997h found that at 1300 K AI 2 O 3 grains are crystalline, while 
iLevin et aP (Il998h reported that the crystalline component only 
dominates at higher temperatures (~ 1450 K). The optical data at 
short wavelengths for AI 2 O 3 are quite uncertain, which translates 
into an uncertainty in the calculated temperature of the amor¬ 
phous AI 2 O 3 grains. For instance, a lower opacity in the near- 
infrared would cause the grains to be cooler. If we consider the 
size and location of the amorphous AI 2 O 3 shell and maintain 
the amounts of grains, the emission produced is proportional 
to the Planck function, since the emitting medium is optically 
thin. Therefore, a lower temperature of 1400 or 1300 K for the 
large grains would correspond to 10 or 20 % less emission, re¬ 
spectively, and the amorphous AI 2 O 3 from the shell could still 
dominate the emission at 12 gm. 

It is also possible, and quite likely, that crystalline and amor¬ 
phous AI 2 O 3 co-exist in the envelope of oxygen-ri ch AGB stars. 
In fac t, grains consisting of pure crystalline AI 2 O 3 (IZeidler et alJ 
l2013h and th ose consisting of a crystalli ne AI 2 O 3 core and a sil¬ 
icate mantle (iKozasa & Sogawall997ah have been proposed as 
the carriers of the unidentified 13 pm feature. We tried to fit the 
13 pm feature of WHya with pure crystalline AI 2 O 3 grains in 
the GB DS using the optical constant measured bv IZeidler et alJ 
(l2013h . We can only obtain a good fit by considering spherical 
(rather than DHS) particles and the optical constants for crys¬ 
talline AI 2 O 3 grains at 973 K. Our models show that the mass 
abundance of crystalline AI 2 O 3 grains compared to the amor- 
pho us ones would be very small, around 5% by mass. We note 
that iTakigawa et al.l (l2014l) found a high fraction of crystalline 
pre-solar AI 2 O 3 grains. The authors measured five out of nine 
studied pre-solar AI 2 O 3 grains, which were probably produced 
in the ejecta of a low-mass AGB star. The derived degree of crys¬ 
tallinity is considerably higher than that inferred from the model 
presented in this work. 

4.4.2. Silicate grains 

We can speculate on what would be the lattice structure of sil¬ 
icate solids if they condensed on a warm AI 2 O 3 core in the 
GBDS. If a silicate solid is heated to a temperature above 
its glass temperature (~ 1000K), it will crystallize. Large pure 
Mg 2 Si 04 grains have a temperature of 700 K at 2.0 R*; a popula¬ 


tion of small pure amorphous AI 2 O 3 grains at the same distance 
is heated to about 1300K. The temperature of a composite grain 
will obviously depend on the relative fraction of each species. If 
a low mass of Mg 2 Si 04 condenses onto the surface of an amor¬ 
phous AI 2 O 3 core the Mg 2 Si 04 solid is expected to be crystalline 
because the temperature of the composite would rise above the 
glass temperature for silicates. However, if 35% silicon conden¬ 
sation and full aluminum condensation is assumed, the silicate 
grain mass would be approximately 11.5 times that of amor¬ 
phous AI 2 O 3 and the composite grains would probably remain 
cold and undetectable. A rough calculation suggests that when 
the mass of Mg 2 Si 04 and amorphous AI 2 O 3 are comparable, the 
temperature of the composite grain would be about lOOOK at 
2.0 R*, that is close to the glass temperature for silicates. This 
means that only roughly the first 10 % of silicate material to con¬ 
dense onto an amorphous AI 2 O 3 grain in the GBDS would yield 
a crystalline solid, while continued silicate condensation would 
make the grain amorphous. 

4.5. Expansion veiocity prof He of the wind of W Hya 

In our best-fit model, the silicate emission is produced in an en¬ 
velope with inner radius of 50 AU. In order for these grains to 
be visible, large amounts of metallic iron must be placed in ther¬ 
mal contact with the silicate grains. We can connect this find¬ 
ing with the gas model from Paper I. In that work, the wind of 
WHya was found to slowly accelerate up to about 5.5 km/s and 
then to quickly accelerate to the final expansion velocity of 7.5 
km/s. Such an additional acceleration requires extra opacity to 
add momentum to the (supersonic) flow. In a dust-driven out¬ 
flow, changes in the dust properties could in principle provide 
such additional momentum. We speculate that the late conden¬ 
sation of metallic iron on the silicate grains might provide the 
extra opacity needed for an increase in the wind speed to its final 
value. 

This interpretation is based on a steady-state picture for the 
inner wind (r < 1000 AU). An alternative explanation is that 
the conditions at the wind-acceleration region have changed 
about 30 years ago, causing the dust composition to change 
from silicate-rich to silicate-poor and the wind expansion ve¬ 
locity to decrease. However, we do not see signs of significant 
mass-loss rate variations in the acceleration zone, as one might 
expect if the local conditions had significantly changed. In ad¬ 
dition, the expansion velocity measured using different diag- 
nostics has been shown to increase stead ily outwards (see, e.g., 
iMuller et al.ll2008[ ISzvmczak et al.lll998h . Moreover, if the ma¬ 
terial at r > 50 AU and now expanding at 7.5 km s ' had been 
ejected with a higher velocity from the start than material in¬ 
wards, one would expect to see a gap (of about 15 AU) between 
these two wind components. There is no support for such a gap 
in the data at hand. Although we are not able to rule out this 
alternative explanation completely, we consider it unlikely. 

4.6. Recent mass-toss history of WHya 

We And that W Hya has undergone abrupt changes of more than 
an order of mag nitude in i ts dus t mass-loss rate in the past few 
thousand years. iHawkinsI (Il99(ll) reported observations carried 
out with IRAS that reveal a dust shell of 1 x 10“^ Mq extending to 
about 15 '. If the same expansion velocity is considered, it would 
take the dust grains forty-five thousand years to reach such dis¬ 
tances from the star and, therefore, the average dust mass-loss 
rate in that period would be 2 x lO^^Moyear *. This value is 
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compatible with the mass-loss rate derived by us for the phases 
with high mass-loss rates (phases II, IV and V, given in Table 
S, rather than to the more recent period of low mass-loss rate 
(phases I and II, given in Table |2l). Therefore, if the mass loss 
fluctuates between a low and high mode, the high mode must 
be the norm. W Hya might have recently undergone a change in 
mass-loss rate that had not occurred in the past fifty thousand 
years. 

Although thin shells produced by similar mass-loss rate dis¬ 
continuities are found somewhat frequently around carbon-rich 
AGB stars (e.g. lOlofsson et alj[l9^ 1199^ . the narrow mass- 
loss rate peak seen between 5400 and 5900 AU in the PACS im¬ 
ages o f W Hya is unusual for oxygen-rich AGB stars (ICox et alJ 
I 2 OI 2 I) . The shells around carbon-stars are thought to form as a re- 
sult of a strong increase in mass-loss rate duri ng a thermal-pulse 
(ISchroederet al.lfl998t iMattsson et ^l2007h . It is not known 
why these detached shells are not observed around oxygen- 
rich AGB stars, and W Hya is unique in that sense. The zero- 
age-mai n-sequence mass of W Hya was determined to be 1.5 + 
0 . 2 Mo (iKhouri et alJl2014H) . which plac es this object around 
the li mit for becoming a carbon-rich star (IWallerstein & Knannl 
Il998h . However, models for the formation of detached shells re¬ 
quire the i nteraction of a faste r wind with a previously slower 
wind (e.g. ISchoier et al.ll 200 ^ . and that is not expected to be 
directly connected to stellar mass, but rather to changes in the 
wind properties during the thermal pulse. Hence, it is currently 
not clear whether the peak in the mass-loss rate of W Hya is due 
to a thermal pulse or whether the star can be seen as an oxygen- 
rich analogue to carbon stars with detached shells. 


5. Summary 

We presented a model for the dust envelope of W Hya that simul¬ 
taneously considers different observ ables of the dust en velope of 
this star: t he scattered light fracti ons jNorris e t alJ2012h . the ISO 
spectrum (lJusttanont et al.ll2004l : ISloa n et al.ll2 003h. th e inner ra¬ 
dius of the silicate emission shell ( Zhao-Ge isler et ani2011|). the 
elemental abundances in the wind (Paper I. iKhouri et al.l2014bh . 
and the PACS 70 pm image (ICox et al.ll2012[) . 

Our model consists of a gravitationally bound dust shell 
combined with an outflow. An MRN distribution of amorphous 
AI 2 O 3 particles in the GBDS fits both the infrared excess emis¬ 
sion and the scattered light fractions. However, large translu¬ 
cent silicates that show no spectral signature at 9.7 pm, but that 
are able to reproduced the scattered light fractions, can exist 
in the GBDS and wind, implying that we are unable to de¬ 
termine the nature of the scattering agents. Silicate grains are 
preferred, however, based on wi nd-driving models because of 
the lo w abundance of aluminum (lH6fneill2008HBladh & Hofneil 
l2012h . Independent of the nature of the grains that drive the 
outflow, spectral emission from the silicate grains arises from 
distances beyond 50 AU from the star. This emission may arise 
when iron inclusions are incorporated into the silicates heating 
up the grains. The silicon content that is in the dust is found to be 
consistent with the amount of silicon atoms that dis appears from 
the gas phase, as extracted from SiO line emission (IKhouri et al.l 
l2014bl) . We find evidence that calcium is a constituent of the dust 
in the wind of W Hya. It is difficult, however, to quantify the cal¬ 
cium content of the dust because we only have optical constants 
available for three calcium-bearing species. Although our best fit 
seems to suggest a super-solar Ca abundance in the outflow, op¬ 
tical data over a broad wavelength range of silicates with lower 
calcium content might result in fits with a solar calcium abun¬ 
dance. Emission from calcium-bearing dust species might also 


originate from the GBDS, as long as there is no strong contribu¬ 
tion in the MIDI wavelength range. 

The mass-loss rate of W Hya seems to have been at the 
present level for about 300 years. At times before this, we see a 
phase when the mass-loss was lower by a factor of roughly two, 
which is traced by both the dust envelope and the CO gas emis¬ 
sion. The mass-loss rate remained at this level for about 2500 
years and produced dust now located in the envelope between 
500 and 4000 AU. This phase of lower mass loss was preceded 
by one where W Hya lost mass at a rate at least three times higher 
than at present, with a localised peak in which the rate was even 
twenty times higher and which ended about 3400 years ago and 
lasted for about 300 years. 
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Table A.l. Diagnostic parameters of the silicate dust species tested by us in the fitting procedure. and 

represent the peak value of the opacity in the 9.7/tm and 18 /tm silicate peaks and the ratio between these two values, respectively. 
A and A are the wavelengths where the two peaks are found. 


Species 

Ref. 

^max 
^9.7 um 

[cmVg] 

1 max 
^ ^9.7 

j^ax 

^18um 

[cmVg] 

A fc-max 

[jum] 

„max /„max 
^18/2m/^9.7/im 

CajAbSiO, 

1 

3210 

10.4 

1270 

20.7 

0.40 

Ca2Mgo,5Al2Sii,507 

1 

3920 

10.3 

1450 

19.4 

0.37 

Mgo,4Feo.6Si03 

2 

3050 

9.5 

1680 

20.9 

0.55 

Mg0.5F6o,43Cao.03Alo,04Si03 

3 

3640 

9.7 

1670 

18.9 

0.46 

Mgo,5Feo,5Si03 

2 

3370 

9.5 

1720 

19.3 

0.51 

Mgo,7Feo,3Si03 

2 

3940 

9.6 

1970 

18.6 

0.50 

Mg0.7SiO2.7 

4 

2720 

9.3 

1530 

20.9 

0.56 

Mgo.8Feo2Si03 

2 

4240 

9.7 

2000 

18.6 

0.47 

Mgo,8Fei,2Si04 

2 

2610 

10.2 

1800 

18.0 

0.69 

Mg,,5Fei,5AlSi30,o.5 

1 

3790 

10.2 

1700 

18.8 

0.45 

Mg|.5Si03.5 

4 

2460 

9.8 

1440 

18.9 

0.59 

Mg2.4Si04.4 

4 

1590 

10.2 

1370 

17.8 

0.86 

Mg2AlSi207,5 

1 

4480 

10.0 

1740 

18.7 

0.39 

Mg2Fe2AlSi40i3.5 

1 

3850 

10.1 

1730 

18.7 

0.45 

Mg2Si04 

4 

2440 

9.9 

1220 

17.7 

0.50 

Mg3AlSi30io.5 

1 

4390 

10.0 

1790 

18.5 

0.41 

Mg4.5 Fe4 5 AlSi9028.5 

1 

3380 

10.1 

1820 

18.5 

0.54 

Mg4AlSi40l3.5 

1 

4450 

9.9 

1810 

18.4 

0.41 

MggAlSigObg.S 

1 

4750 

9.8 

1980 

18.3 

0.42 

MgFeAlSi267.5 

1 

3640 

10.3 

1590 

18.6 

0.44 

MgFeSi04 

2 

2700 

10.1 

1873 

18.0 

0.69 

MgSi03 

2 

4910 

9.6 

2180 

18.4 

0.44 


Notes. 1 - iMutschke et all (Il998li . 2 - iDorschner et all (Il99.'5h . 3 - Llageret^ 119941) . 4 - Llager et ^ l200.3f) 
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